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ABSTRACT 

Using the standard dynamical theory of spherical systems, we calculate the properties 
of spherical galaxies and clusters whose density profiles obey the universal form first 
obtained in high resolution cosmological TV-body simulations by Navarro, Frenk & 
White. We adopt three models for the internal kinematics: isotropic velocities, con- 
stant anisotropy and increasingly radial Osipkov-Merritt anisotropy. Analytical solu- 
tions are found for the radial dependence of the mass, gravitational potential, velocity 
dispersion, energy and virial ratio and we test their variability with the concentration 
parameter describing the density profile and amount of velocity anisotropy. We also 
compute structural parameters, such as half-mass radius, effective radius and various 
measures of concentration. Finally, we derive projected quantities, the surface mass 
density and line-of-sight as well as aperture velocity dispersion, all of which can be 
directly applied in observational tests of current scenarios of structure formation. 

On the mass scales of galaxies, if constant mass-to-light is assumed, the NFW 
surface density profile is found to fit well Hubble-Reynolds laws. It is also well fitted 
by Sersic R 1 /" 1 laws, for m ~ 3, but in a much narrower range of m and with much 
larger effective radii than are observed. Assuming in turn reasonable values of the 
effective radius, the mass density profiles imply a mass-to-light ratio that increases 
outwards at all radii. 



Key words: 

of Universe 



methods: analytical - galaxies: clusters: general - large-scale structure 



1 INTRODUCTION 

A universal profile of dark matter haloes was introduced as 
a result of high-resolution A^-body simulations performed 
by Navarro, Frenk & White (1995, 1996, 1997, hereafter 
NFW) for power-law as well as CDM initial power spectra 
of density fluctuations. NFW found that in a large range 
of masses the density profiles of dark haloes can be fitted 
with a simple formula with only one fitting parameter. The 
density profile steepens from r _1 near the centre of the halo 
to r -3 at large distances. The NFW profile has been con- 
firmed in cosmological simulations by Cole & Lacey (1996), 
Tormen, Bouchet & White (1997), Huss, Jain & Steinmetz 
(1999a), Jing (2000), Bullock et al. (1999), while Huss, Jain 
& Steinmetz (1999b) have shown that the NFW profile also 
arises from non-cosmological initial conditions. It is worth- 
while noting that some (but not all) recent very high resolu- 
tion cosmological simulations produce steeper density pro- 
files, with inner slopes ~ —1.5 (Fukushige & Makino 1997, 
Moore et al. 1998, Ghigna et al. 1999, see also Jing & Suto 
2000). The density profiles in the cosmological simulations 



also display considerable scatter (Avila-Reese et al. 1999, 
Bullock et al. 1999), and Avila-Reese et al. find that the 
outer slopes of galaxy size haloes are steeper than the NFW 
slope of —3 when selected within clusters (—4) and slightly 
shallower within groups (—2.7). Although the exact proper- 
ties of dark matter haloes are still under debate, the NFW 
profile is presently considered to provide the reference frame 
for any further numerical research on density profiles of dark 
haloes. 

Simple cosmological derivations of the density profiles 
of bound objects are difficult, essentially because one needs 
to work in the non-linear regime of the growth of gravi- 
tational instabilities. Nevertheless, using the spherical top- 
hat model of Gunn & Gott (1972), density profiles typically 
varying as r _9//4 were derived by Gott (1975), Gunn (1977), 
Fillmore & Goldreich (1984) and Bertschinger (1985). Hoff- 
man & Shaham (1985) applied the spherical infall model to 
the hierarchical clustering scenario and predicted that the 
density profiles of haloes should depend on Q as well as the 
initial power spectrum of density fluctuations. However, for 
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fl = 1 they obtained power-law profiles in contradiction with 
the steepening slopes found in the current TV-body simula- 
tions described above. In a recent study, Lokas (2000) has 
improved the model of Hoffman & Shaham (1985) by a gen- 
eralization of the initial density distribution, the introduc- 
tion of a cut-off in this distribution at half the inter-peak 
separation and by a proper calculation of the collapse fac- 
tor. The improved model reproduces the changing slope of 
the density profile and its dependence on halo mass and the 
type of cosmological power spectrum found by NFW. The 
NFW profile is also reproduced in studies taking into ac- 
count the merging mechanism (see Lacey & Cole 1993) in 
the halo formation scenario (e.g. Salvador-Sole, Solanes & 
Manrique 1998, Avila- Reese, Firmani & Hernandez 1998). 
Therefore the numerical and analytical considerations seem 
to converge on the statement that the density profiles of 
dark matter haloes are indeed well described by the univer- 
sal formula proposed by NFW. 

The ultimate test of both the analytical and numerical 
results must come from the observations of density profiles 
of galaxies and galaxy clusters. Three recent studies of clus- 
ters (Carlberg et al. 1997, Adami et al. 1998, van der Marel 
et al. 2000) claim good agreement between cluster obser- 
vations and the NFW mass density profile. But for galax- 
ies, the situation is less satisfying. Flores & Primack (1994) 
show that the NFW profile is incompatible with the rota- 
tion curves of spiral galaxies, while Kravtsov et al. (1998) 
estimate that the inner slope of the density profile of dwarf 
irregular and LSB galaxies is —0.3 instead of —1. However, 
these conclusions were obtained with a number of assump- 
tions and approximations concerning the very unclear issues 
of biasing, non-sphericity of objects and so on. Besides, as 
pointed out by van den Bosch et al. (2000), Swaters, Madore 
& Trewhella (2000) and van den Bosch & Swaters (2000), 
the observed rotation curves of these galaxies are too uncer- 
tain to discriminate between cores and cusps. 

The main motivation for this research is to explore ana- 
lytically the physical properties of objects with NFW density 
profiles. The aim is to check whether these properties are ac- 
ceptable from the physical point of view and thus to test the 
validity of density profiles obtained in cosmological TV-body 
simulations. Additionally, this paper presents formulae for 
observable quantities that can be used for comparisons be- 
tween the theoretical predictions (such as the NFW profile) 
and observations. 

The paper is organized as follows: after a short presen- 
tation of the universal formula for the density profile pro- 
posed by NFW, in Section 2 we describe physical proper- 
ties of spherical systems following from this density profile. 
Section 3 is devoted to a simple comparison between the 
projected NFW density profile and the surface brightness of 
elliptical galaxies. A more thorough comparison is beyond 
the scope of the present paper and will be given elsewhere 
(Mamon & Lokas, in preparation). The discussion follows in 
Section 4. 



2 PROPERTIES OF THE NFW MODEL 
2.1 Basic properties 

NFW established that the density profiles of dark matter 
haloes in high resolution cosmological simulations for a wide 



range of masses and for different initial power spectra of 
density fluctuations are well fitted by the formula 



p(r) 



(r/r s ) (1 + r/r s 



(1) 



with a single fitting parameter <5 c har, the characteristic den- 
sity. The so-called scale radius r s is defined by 

= (2) 

where r v is the virial radius usually defined as the distance 
from the centre of the halo within which the mean density 
is v times the present critical density, p°. The value of the 
virial overdensity v is often assumed to be v = 178, a number 
predicted by the simplest version of the spherical model for 
= 1. For other cosmological models it can be lower by a 
factor of 2 or more (Lacey & Cole 1993, Eke, Cole & Frenk 
1996). However, according to the improved spherical infall 
model (Lokas 2000) v can be as low as 30 even for Q = 1. 
In the following, v is kept as a free parameter. 

The quantity c introduced in equation (^) is the con- 
centration parameter, which is related to the characteristic 
density by 



vc 3 g(c) 



where 



9(c) 



ln(l + c) - c/(l + c) 



(3) 



(4) 



The concentration parameter will be used hereafter as the 
only parameter describing the shape of density profile. From 
cosmological TV-body simulations (Navarro et al. 1997, Jing 
2000, Bullock et al. 1999, Jing & Suto 2000), extended Press- 
Schechter theory (Navarro et al. 1997, Salvador- Sole, Solanes 
& Manrique 1998), and the spherical infall model (Lokas 
2000), we know that c depends on the mass of object and the 
form of the initial power spectrum of density fluctuations. 
For all initial power spectra, the observed trend is for lower 
concentration parameter in higher mass objects, with 4 < 
c < 22 in cosmological simulations with CDM initial power 
spectra and c up to 90 for the less realistic scale-free power 
spectra. More precisely, in the ACDM cosmology, c = 5 
corresponds to the masses of clusters of galaxies, while c = 
10 corresponds to the masses of bright galaxies. 

It is convenient to express the distance from the centre 
of the object in units of the virial radius r v : 



and the density profile of equation (jjj) then becomes 

P(s) = vc 2 g(c) 
p° 3s(l + cs) 2 ' 



(5) 



(6) 



The mass of the halo is usually defined as the mass 
within the virial radius: 

M v = -nrlvp° . (7) 

The distribution of mass in units of the virial mass follows 
from equation (ph: 



M(s) 
M v 



= g(c) 



ln(l + cs) 



(8) 
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and we see that it diverges at large s, which is a disadvantage 
of the model from a physical point of view. 

The gravitational potential associated with the density 
distribution Q is 



*(«) 

Vv 



-3(c) 



ln(l + cs) 



where V v is the circular velocity at r 
GM V 4 



ri 2 
7r G r„ v p c 



(9) 



(10) 



Hence, from equation (b|) we see that the gravitational po- 



tential at the centre, $(0 
Equations (^) and ([" 
obeys 

F 2 (s) 



9(c) 



V 2 



ln(l + cs) 
(0) and (T 



■ —cg(c)V 2 , is finite, 
lead to a circular velocity that 



(11) 



1 + cs. 

) were first derived by Cole & 



Equations 
Lacey (1996). ~ 

The radial velocity dispersion a T (r) can be obtained by 
solving the Jeans equation 

2 



d£ 
dr 



(12) 



p dr 

where (3 = 1 — a 2 (r)/a 2 (r) is a measure of the anisotropy 
in the velocity distribution. In the simplest case of isotropic 
orbits, ag(r) = cr r (r) and (3 = 0. This value of /3 is also 
close to the results of Af-body simulations: Cole & Lacey 
(1996) and Thomas et al. (1998) show that, in a variety of 
cosmological models, the ratio /Je/tJ r is not far from unity 
and decreases slowly with distance from the centre to reach 
~ 0.8 at the virial radius. However, Huss, Jain & Steinmetz 
(1999a) find (Je/o, ~ 0.6 at r v . 

First we consider the case of /3=const. Then the solution 
to the equation ( |l2| ) with the condition of o l . — *• at s — *• oo 
is 

.2 

const) = g(c)(l + cs) 2 s 1-2 ' 3 



th( s ^ 



ln(l + cs) cs 



2/3-2 



(1 +cs) 2 



(1 + cs) 



ds. (13) 



For f5 = 0, 0.5 and 1, reasonably simple analytical solutions 
to this equation can be found: 

2 



a 

v;? 



1 



\c g(c) s (1 + cs) 2 [tt 2 



ln(cs) 

cs 



(1 + cs) 



1 + cs 



+ 1 



1 



cs 1 + cs 



x ln(l + cs) + 3 In (1 + cs) + 6Li 2 (-cs) 



(14) 



C )(l + cs) 2 [- 



+ 



2(1 + ca) 



+ 



ln(l + cs) + ln(l + cs) 



1 + cs cs 1 + cs 

- ln 2 (l + cs) - 2Li 2 (-cs)], 



(15) 



In the above expressions Li2(a;) is the dilogarithm, a 
special function which can be conveniently dealt with using 
Mathematica packages. Otherwise, it can be approximated 
by 



Lh(x) = 



ln(l-t)dt 



t 



'x[l+lQ-°- 6 (-x) 



0.62/0.7] 



.(17) 



The fit is accurate to better than 1.5% in the range —100 < 
x < 0. 

We included the predictions for j3 — 1 just as a limiting 
case. In fact such a model with purely radial orbits and 
NFW density profile is not physical since its distribution 
function is not everywhere non-negative. As pointed out by 
e.g. Richstone & Tremaine (1984, see also Lokas & Hoffman 
2000), such velocity anisotropy requires the inner density 
profile to be r~ 2 or steeper for the model to be physical. 

A more realistic description of velocity anisotropy is 
provided by a model proposed by Osipkov (1979) and Mer- 
ritt (1985) with f3 dependent on distance from the centre of 
the object 



Pom 



s 2 + si 



(18) 



where s a is the anisotropy radius determining the transi- 
tion from isotropic orbits inside to radial orbits outside. As 
mentioned above, the results of A^-body simulations suggest 
a$/(T r — 0.8 and therefore /3 ~ 0.36 at s = 1, which gives 
s a ~ 4/3, a value that we adopt here for all numerical cal- 
culations. 

For the Osipkov-Merritt model the solution to the Jeans 
equation (with the same boundary condition as before) reads 



t^(s,/3om) 

Vv 



g(c)s(l + cs) 2 
s 2 + s 2 

(s 2 + s 2 )ln(l + cs) c(s 2 + s 2 ) 



'(l + cs) 



and the integration gives 



s 2 (l + cs) 3 



ds (19) 



t^-(s,/3om) 



CSa 

s 



g (c)s(l + cs) 2 
2(s 2 + sl) 

c 2 s\ ln(cs) + c 2 s 2 ln(l + cs) ( 1 - 



(l + c 2 s 2 ) 



+ (l+3c 2 s 2 



1 



+ 



2 ln(l + cs) 



(1 + cs) 2 

2__ 

3 1 + cs 



(20) 



■ In 2 (1 +cs) +21+2 (-cs) 



Figure ^ shows the radial dependence of the radial ve- 
locity dispersion. The upper panel of the Figure presents 
how the results depend on the concentration parameter in 
the isotropic case, while the lower panel compares predic- 
tions for different anisotropy models with c = 10. 



1 



9(c) (1 + cs) 2 



2(1 + cs) 2 



1 + cs 



Hl + cs) ln 2 (l + cs) ' 
-i , — H 7, hLi 2 (-cs) 



.(16) 



© 0000 RAS, MNRAS 000, 000-000 



4 Ewa L. Lokas and Gary A. Mamon 




_i i i i i i 



0.01 



0.01 





r/r v 



r/r v 



Figure 1. Radial velocity dispersion profile (in units of the cir- 
cular velocity at the virial radius), given by in isotropic model, 
equation (|14[) , for three different values of the concentration pa- 
rameter c (upper panel) and for the four considered anisotropy 
models with c = 10 (lower panel). 



Figure 2. The radial dependence of the virial ratio in the 
isotropic model (eqs. plj and j^J) for three different values of 
the concentration parameter (upper panel) and for the four con- 
sidered anisotropy models with c = 10 (lower panel). 



2.2 The energy distributions 

The potential energy associated with the mass distribution 
of equation (H) is 



W(s) = 

r v 

= -Wo, 

where 

W aa = - lim W(s) 



GM(s) dM(s) 



ds 



s ds 

1 21n(l + cs) 



(1 + cs) 2 
cg 2 (c)GM 2 



1 + cs 



2r v 



The kinetic energy for arbitrary ft is given by 

T(s,P) =2-ktI [ (3-2/3)p(s)a^(s,f3)s 2 ds. 
Jo 



(21) 



(22) 



(23) 



For the three cases of 0=0, 0.5 and 1, we obtain respectively 

T{s,f3 = 0) = ivy oo {-3 + T -^-21n(l + cs) 

+ cs[5 + 3 ln(l + cs)] - cs 2 [7 + 6 ln(l + cs)] 
+ c 3 s 3 [7r 2 — In c — In s + ln(l + cs) 



+ 3m 2 (l + cs) + 6Li 2 (-cs)]} 



T(s,/3 = 0.5) = i^ 0O {-3 + T -^-31n(l + cs ) 
+ 6cs[l + ln(l + cs)] — c 2 s 2 [7r 2 



+ 31n 2 (l + cs)+6Li 2 (-cs)]} 



(24) 



(25) 



T{8,0=1) = iw 00 {-21n(l + cs) + cs[^-- T -J— 
I 6 1 + cs 



+ ln 2 (l + cs) +2Li 2 (-cs)]}, 



(26) 
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Figure 3. Dependence on the concentration parameter of the 
virial ratio at the virial radius for the four considered anisotropy 
models. 



Figure 4. Dependence on the concentration parameter of the 
half-mass radius, scaled to the virial radius (thicker lower solid 
line, see eq. [Q) and 7 for the four anisotropy models (eq. Jjo| , 
four upper lines). 



where we have used in each case the corresponding expres- 
sion for a 2 (s,f3) from equations (|L4^-(|l6|). For the Osipkov- 
Merritt model the calculation has to be done numerically. 

The results for the potential and kinetic energy (|2l|)- 
( p6[ ) lead to a virial ratio lim,,^^ 2T/\W\ = 1 for any value 
of c, in agreement with the virial theorem. Figure pj shows 
how the virial ratio depends on distance for three different 
values of the concentration parameter in the isotropic case 
(upper panel) and compares the ratios obtained for different 
(3 with c = 10 (lower panel). At low radii, the virial ratio is 
large, especially for low concentration parameters and mod- 
els with much anisotropy. However, as demonstrated by Fig- 
ure H at the virial radius r„(s = 1), 2T/| W\ is still greater 
than unity and grows with the amount of anisotropy in the 
model. We see that the virial theorem is better satisfied at 
s = 1 for objects with larger concentration parameters, as 
lim c ^oo 2T/\W I (s = 1) = 1. Since objects of smaller mass 
have larger concentration parameters, they are closer to dy- 
namical equilibrium. 

The scalar virial theorem we referred to above is ex- 
pected to be satisfied for self-gravitating systems in steady 
state. In more realistic situations, the system is never iso- 
lated and experiences an external gravitational field; there is 
also continuous infall of matter. We may conclude from the 
results above that objects with NFW density profiles and 
different velocity distributions are close to dynamical equi- 
librium. However, the virial ratio cannot be used to define 
the boundary of the virialized object. 



2.3 Structural parameters 

A useful quantity is the half-mass radius. Unfortunately, the 
divergence of the mass of the NFW profile forces one to de- 
fine the half-mass radius within a cutoff radius r cut . The 
most natural choice is r cut = r v , since the density distribu- 
tion is only reliable out to the virial radius. With r cut = r v , 



the half-mass radius satisfies the following relation for the 
mass of dimensionless radius: 



M(l) 



(27) 



Numerical values of rh/r v are easily obtained using equation 
(^) and over the range 1 < c < 100 they can be approxi- 
mated to better than 2% accuracy by 



0.6082 - 0.1843 logc 
- 0.1011 log 2 c + 0.03918 log 3 c 



(28) 



The lowest thick solid line in Figure |j shows how rh/r„ 
decreases with increasing concentration parameter. 

It is useful to estimate the concentration 7 of a dynam- 
ical system, such that 



/ 2\ GM 
(a \ = 7 

v 1 r h 



(29) 



where {a 2 } = (a 2 + u'g + a'^} is the mass weighted mean- 
square velocity dispersion. As first noted by Spitzer (1969) 
for polytropes, many realistic density profiles have 7 = 0.4. 
For example, it is easy to show that for the Hernquist (1990) 
model with (3 = 0, 7 =<1 + y/2)/6 ~ 0.403 (Mamon 2000). 

Using equation (E9I) and limiting again the mass to 
r C ut — r v , we define 7 with 



Hi 



<* 2 > 



_ _ r h T(l,/?) 

' GM(1) * GM 2 (1) 

where we made use of 
T[x,P) = \M{x) <- 2 ) r< _ • 



(30) 



(31) 



The values of 7 for different velocity anisotropy models, de- 
rived from equations (0), (§), (||), (§§)-©, (f|), and (|o|) 
are shown in Figure k and in the case of f3 = yield num- 
bers closest to 0.4: 7 = 0.56 for c — 5 and 7 = 0.51 for 
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c = 10. Thus the NFW model produces 7s that are higher 
than the canonical value of 0.4, especially if more velocity 
anisotropy is assumed. This may be caused by the ill-defined 
cutoff radius. 

In models with homogeneous cores, the central density, 
the core radius r c and the central 3-D velocity dispersion 
a 2 (0) are related through 



(32) 



King (1966) models have 77 — 9. In models with cuspy cores, 
we propose the scaling relation 



4nGp(r s )r 2 = \v(° 2 ) r 



Using equations (121) 



(33) 



and 



one has 4nGp(r s 



cg(c) V 2 /4 and from equation (fyh for x = 1/c one obtains 
3cg(c)V 2 M(l/c) 



V 8T(l/e,/?) • 
For different velocity anisotropy models we then have 
3(2 In 2 - 1 



V(P = 0) 



VW = 0.5) 



1,09 = 1) = 



2(tt 2 -7-81n2 + 61n 2 2) 

9(1 -2 In 2) 
~ 4(tt 2 - 9-61n2 + 61n 2 2) 

9(21n2-l) 



2(tt 2 



12 In 2 + 6 In 2 2) 



2.797, 



2.138, 



1.212, 



(34) 



(35) 



(36) 



(37) 



where we have used equations (g) and (|24j)-(]26]), and the fact 
that Lia(— 1) = — 7r 2 /12. Note that 77 is independent of c in 
all cases with (3 =const. For the Osipkov-Merritt model rj is 
no longer a constant but we find 1.902 < 77 < 2.797 in the 
range 1 < c < 100 with the limiting cases of 77 — > r)(/3 = 1) 
for c — > and 77 — > 77(73 = 0) for c — *• 00. Such limiting 
behaviour is due to the fact that for large c the integration 
of T(l/c, /3), equation (p3~|), probes only the range of s where 
/3 is close to zero, while for small c the integral is dominated 
by contribution from large s where /3 is close to unity. 
Finally, we consider the structural parameter 



WUM : 



W(s) 
M(s)<&(0) 



(38) 



brought forward by Seidov & Skvirsky (2000) with the moti- 
vation of WUM being constant for different self-gravita ting 
objects of simple geometry. Using equations (^) , and (J2TJ) 
we find that for the NFW model 



WUM : 



cs(2 + cs) - 2(1 + cs) ln(l + cs) 
2(1 + cs)[-cs + (1 + cs) ln(l + cs)] 



(39) 



so the parameter turns out to be a function of cs — r /r s only. 
It grows with s from zero at s — ► reaching a maximum 
value of 0.196 at r/r s = 4.62 and decreases to zero again as 
s — > 00. The values of this parameter at the virial radius are 
0.196, 0.187 and 0.125 respectively for c = 5, 10 and 100. 

2.4 The distribution function 

A quantity of great dynamical importance is the distribu- 
tion function. For a spherical system with an isotropic veloc- 
ity tensor, the distribution function depends on the phase- 
space coordinates only through the energy (e.g. Binney & 




_ D I 1 1 1 1 I 1 1 1 1 L_ 

-12.5 -2.0 -1. 



_i L_i 1 1 1 I 1 ' ' 

5 -1.0 -0.5 0.0 



Figure 5. The distribution function for isotropic model (eq. |40| ) 
for three different values of the concentration parameter. 



Tremaine 1987), and can be derived through the Eddington 
(1916) formula (e.g. Binney & Tremaine 1987): 



V8tt 2 



A 2 p d* 1 /dp\ 

d* 2 VF^* + Vdl>J*=o 



(40) 



where £ and \P are the conventionally defined relative energy 
and potential; here £ = — E, where E is the total energyjaer 
unit mass and 9 = — $, where $ is given by equation 

It is easy to show that, given equations (jfjj) and ^ 
the second term in brackets in equation (^) is zero. The 
simplest way to perform the integration of the first term is to 
introduce dimensionless variables ^ = VP/Ci and p = pjCi, 
where Cj = g(c) V 2 and C 2 = c 2 g(c)M v / (Anri). Then the 
integration variable should be changed to s and the limit of 
integration corresponding to £ found numerically for each £ 
by solving equation \f/ (s) = £ . Otherwise, with a few percent 
accuracy, the integration in (^o|) can be done directly with 
an approximation s apx = — 1. 751n(\I//c) /W. 

The calculations of the distribution function are usually 
performed in units such that G = M — R c — 1 (Binney & 
Tremaine 1987), where M is the total mass of the system and 
R c is its effective radius. Since in the case of NFW profile 
the total mass is infinite a reasonable choice seems to be to 
put M v — 1. The effective radius is not well defined either 
but can be approximated as r v /2 (see the next section). 
Therefore we choose the units so that G — M v = r v /2 = 1 
and arrive at the numerical results shown in Figure Bl This 
choice of normalization is equivalent to measuring / in units 
of VSM v /(r v V v ) 3 and E in units of U„ 2 . 

Figure ^ proves that the distribution function turns out 
to be similar to the distribution functions obtained from 
other density profiles (see e.g. Figure 4-12 in Binney & 
Tremaine 1987), except that the NFW distribution functions 
do not display the cutoff at nearly unbound energies charac- 
teristic of King (1966) models. The results shown in Figure ^ 
indicate a proper behaviour of the distribution function (it 
is nowhere negative). Quantitative comparisons with other 
models should, however, be made with caution because of 



© 0000 RAS, MNRAS 000, 000-000 



Properties of galaxies and clusters 7 




0.0 1 — — — 

0.01 0.1 1 10 

R/r v 

Figure 6. Upper panel: radial dependence of the line-of-sight ve- 
locity dispersion for isotropic orbits (eq. |4i[ ) on the projected ra- 
dius for three values of the concentration parameter. Lower panel: 
comparison of the line-of-sight velocity dispersion profiles for four 
anisotropy models calculated with c = 10. 
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Figure 7. Upper panel: radial profiles of the aperture velocity 
dispersion in the isotropic model for three concentration parame- 
ters. Lower panel: comparison of the aperture velocity dispersions 
for four anisotropy models calculated with c = 10. 



the aforementioned problem with normalization. Distribu- 
tion functions for more realistic velocity dispersion models, 
like the Osipkov-Merritt model, were recently considered in 
detail by Widrow (2000). 

2.5 Projected distributions 

Of primary importance for comparisons with observations 
are the projected distributions. The surface mass density of 
an object is obtained by integrating the density along the 
line of sight: 



SmOR) = 2 



r p(r) 



( r 2__R2)i/: 



■ dr 



2tt 



c) M v 1 

2 



\c 2 r 2 - 1| - 1 / 2 c- 1 [i/( c je)] 

{c 2 R? - l) 2 



(41) 



where 



C-\x) = 



cos (x) if R > r s 
cosh -1 (a;) if R < r s 



(42) 



In the above expressions R is the projected radius and 
R = R/r v . For the singular case R = r s the independent ex- 
pression in equation (Mil) equals 1/3 and we have T,m(R) = 
c 2 g(c)M v / (Q-nr 2 ,) . An analytical formula equivalent to equa- 
tion ([II]) was derived independently by Bartelmann (1996). 
The projected mass is then given by 



M P (R) = 2tt / RT, M (R)dR 



9(c) M v 



C-y./(cR)] nn (cR 

\ C 2 R 2 _ ill/2 



(43) 



which is logarithmically divergent at large R. C 1 (x) is 
again given by equation (|42|). 
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Another important projected quantity is the line-of- 
sight velocity dispersion which for a spherical non-rotating 
system is (Binney & Mamon 1982) 



Oios(-R) = 



pa, 



,P)r 



dr 



(44) 



Zm{R) J r \~ " r 2 J Vr 2 - R 2 

where Em(B) is given by equation ( fill ) and the radial ve- 
locity dispersions a r (r,/3) for our four models are given by 



equations i 
one has 



|)-@ and @. For circular orbits, 



£m(#) 



(?) 



P^ 2 



7? 2 



: dr 



0, and 



(45) 



where 1/ is the circular velocity given by equations (^) and 
([n]). The upper panel of Figure]^ shows the profiles of line- 
of-sight velocity dispersion (with isotropic orbits), obtained 
through numerical integration of equation (^) for different 
concentration parameters. The lower panel of Figure ^ com- 
pares the radial profiles of line-of-sight velocity dispersions 
obtained for c = 10 for different velocity anisotropy models. 

For more distant or intrinsically small galaxies, as well 
as for groups and clusters, spectroscopic observations are 
often limited to a single large aperture centred on the object. 
The mean velocity dispersion within an aperture (hereafter, 
aperture velocity dispersion) is 



where 



S 2 (R) 
M P (R) 



S 2 (R) = 2tt / E M (P)<rL(P)PdP. 



(46) 



(47) 



In the above expressions R is the radius of the aperture, 
Em(-P) is the surface mass distribution, equatio n ([4 11) , and 
M P (R) is the projected mass given by equation (|43jK 

Inserting the expression for <ti os (eq. Q|) into equation 
([l7j), we obtain a double integral, which after inversion of 
the order of integration is reduced to an easily computable 
single integral: 



S 2 (R) = cg{c)M v 



+ 



2 (s,(3) (s 2 



10 

£2)1/2 



(1 +cs) 2 



of (s,P) s 
(1 + cs) 2 

' f3(R 2 + 2s 2 
3s 2 



(»-¥ 



ds 



- 1 



ds 



(48) 



where as before, R = R/r v , s = r/r v and af(s, 0) for differ- 
ent (3 are given by equations (p^)-(|l6[) and (po|). An alog ous 
expression for circular orbits can be obtained from (Uq) by 
replacing a 2 by V 2 , keeping only the terms proportional to 
(3 and dividing by (—2/3). 

Figure fj] displays the radial profiles of aperture velocity 
dispersion, computed numerically from equation (^8|). From 
the upper panel of the Figure we see that in the isotropic 
case the dependence of the results on the concentration pa- 
rameter is rather strong and monotonic for a given R. The 
lower panel of the Figure compares the predictions for dif- 
ferent velocity anisotropy models. 



3 COMPARISON WITH OBSERVATIONS 

Comparisons of the surface mass density to surface bright- 
ness observations are usually performed with the assumption 
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Figure 8. Radial profiles of the surface mass density, given by 
equation (fil|), (upper panel) and the projected mass, equation 
(M3h, (lower panel) for three different values of the concentration 
parameter. Hubble-Reynolds fits from equation (jisj) are shown 
as thin curves (R H R/r v = 0.119, 0.0640 and 0.00743 for c = 
5, 10 and 100, respectively). For c = 100, the NFW surface mass 
density is virtually indistinguishable from the best-fitting Hubble- 
Reynolds law. 



of constant mass-to-light ratio T = const. This assumption 
is not likely to be physical, because of the different physics 
involved in the assemblies of the dark matter and baryonic 
components of galaxies. In particular, the baryons in ellip- 
tical galaxies may well settle at an early epoch, within a 
radius that is the lower of the radius with virial overden- 
sity v ~ 200 and the radius at which gas can cool to form 
molecular clouds and later stars. The baryons in ellipticals 
will then sit today in a region of overdensity v ^> 200, and 
one then expects T to rise with r, at least at large radii. 

Nevertheless, for simplicity, we check whether the ob- 
servations of elliptical galaxies are consistent with the idea 
that stars are distributed within elliptical galaxies accord- 
ing to the NFW density profile, characterized by a virial 
radius where the mean overdensity is 200. Such a situation 
may arise if the dark matter were negligible within elliptical 
galaxies or distributed precisely like the luminous matter. 
In a forthcoming paper (Mamon & Lokas, in preparation), 
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Figure 9. The dependence of the effective radius, defined in equa- 
tion |5(i|), on the concentration parameter, with various choices 
of -Rcut ■ 



we will check in more detail whether the observations of el- 
liptical galaxies are compatible or not with NFW density 
profiles for the mass distribution. 

For constant mass-to-light ratio we have Em(-R) = 
TI(R), where I is the surface brightness. The radial pro- 
files of I = Eji/T and M p are shown in Figure |[ Both 
quantities are normalized to their values at the virial ra- 
dius. Figure ^ shows that the surface mass density depends 
weakly on the concentration parameter, especially at larger 
distances from the centre. 

Since the surface mass density (eq. jll]) behaves as 
1/R 2 at large distances, one may therefore compare it with 
the Hubble- Reynolds formula (Reynolds 1913), which was 
the first model used to describe the surface brightness pro- 
files of elliptical galaxies: 



(1 + R/Rhr) 2 



(49) 



Rhr is the characteristic radius of the distribution, where 
the surface brightness falls to one-quarter of its central value. 
The thin curves of Figure |^ show that the surface mass den- 
sity of the NFW model (eq. |4l)) is very well fitted by equa- 
tion (^) and the best-fit values of J?hr = Rrb,/t v are 0.119, 
0.0640 and 0.00743 respectively for c = 5, 10 and 100. 

The surface brightness profiles of astrophysical objects 
are often scaled with the effective radius, which we denote 
R e , where the projected luminosity is half the total luminos- 
ity. Given the divergence of the projected mass, we are forced 
again to introduce a cut-off at some scale -Rcut = Rcut r v . We 
then have 




(i?/J?W4 



Figure 10. Upper panel: surface brightness profiles (eq. |4l[ ) for 
three concentration parameters and R cu t = 1. Lower panel: the 
dependence of the surface brightness profiles on the cut-off R C ut 
for c = 10 and R C ut = 3, 3.5, 4, 4.5, and 5 (bottom to top curves). 
In both panels, the i? 1 / 4 law (eq. |p^|) is shown as long dashed 
lines. The vertical lines represent the virial radius (for the three 
concentration parameters in the upper panel and for the 5 values 
of -R cu t in the lower panel, with R cu t increasing from right to 
left). The circles are the data points for the galaxy NGC 3379. 



Figure g shows the effective radius, calculated numerically 
from equations (|l3|) and ( po] ) . For i? C ut = 1, a useful approx- 
imation, good to better than 2% relative accuracy, is: 



R e /r v = 0.5565 - 0.1941 log c 



0.0756 log 2 c + 0.0331 log 3 c 



(51) 



The prediction for the surface brightness I — £m/Y 
with Em given by equation ([ll]) expressed in terms of the 
effective radius and the corresponding effective brightness 
I c — I(R C ) is shown in the upper panel of Figure [H^ for 
different values of the concentration parameter c. For com- 
parison, we also show the de Vaucouleurs (1948) R 1 '* law 
describing the observed surface brightness distribution in gi- 
ant elliptical galaxies: 



I(R) = I a exp{-b[(R/R c ) 1/4 - 1]} 



(52) 



M P (R C ) = M p (R cut )/2 



(50) 



where b — 7.67. Clearly, the NFW surface brightness profiles 
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Figure 11. Comparisons of c = 10 projected NFW models (using eq. Jilt) to Sersic models (eq. |^). The curves represent the NFW 
models (for equally spaced values of i? C ut within the interval indicated in each plot, with R C ut increasing upwards on the left portion of 
each plot). The Sersic law is shown as long dashed lines. The vertical lines represent the virial radius (with -Rcut increasing from right to 
left). 



are poorly fitted by the R 1 ^ 4 law, when using R cu t = r v to 
define the effective radius of the NFW profile. 

The lower panel of Figure [Hi] shows how^the results de- 
pend on the choice of cut-off for c — 10 and R cu t = 3, 3.5, 4, 
4.5, and 5. At first glance, it seems that the NFW profile is 
well fitted by the R 1 ' 4 law, especially for R cu t — 4. However, 
the range of surface mass densities where the fit is excellent 
is roughly 10 2 , and the fit is adequate for a range smaller 
than 10 3 . In contrast, the surface brightness profile of the 
nearby giant elliptical galaxy NGC 3379 (M 105) follows 
the R 1/4 law in a range of 10 magnitudes (de Vaucouleurs 
& Capaccioli 1979), i.e. a factor 10 4 in intensity. 

In order to see how good is the de Vaucouleurs's fit in 
this case in both panels of Figure llOl we plotted a number 
of data points equally spaced in R 1 ^. Since de Vaucouleurs 
& Capaccioli (1979) do not provide the error bars for their 
data, the error bars shown in the Figure were taken from 
Goudfrooij et al. (1994). The excess of the data above the 
R 1/4 law for small R was already noted by de Vaucouleurs & 
Capaccioli (1979). The error bars are negligible for R < R c 
and smaller than 15% out to 2.5R a , the maximum distance 
from the centre reached in the data of Goudfrooij et al. 
(1994). 



According to de Vaucouleurs & Capaccioli (1979), in 
this galaxy the R 1/l4 surface brightness profile extends to 
Rum = 7.5 R c = 26.4 kpc, given a distance of 12.4 Mpc 
to NGC 3379 (Salaris & Cassisi 1998). Within i? Um , de 
Vaucouleurs & Capaccioli (1979) report a blue magnitude, 
corrected for galactic extinction of B = 10.10, yielding a 
total blue luminosity of 2.2 x W 10 L@, hence a blue lumi- 
nosity density of 2.8 x 10 s Lq kpc -3 . Since the mass within 
Riim must be greater than the mass in stars, we infer that 
within this radius, Tb > 8 (the typical mass to blue lumi- 
nosity ratio for old stellar populations), yielding an over- 
density of the galaxy, relative to the critical density p c 
of v > 1.6 x 10 4 /(-ffo/70kms _1 Mpc" 1 ) 2 . Therefore, since 
v 2> 100 (the value at r v ), we conclude that Rn m <C r v , 
hence i? e <C r v /7.5. In contrast, with -R cu t = 1, the effective 
radius of the NFW model (c = 10) is ~ 0.3 r v (Figure |). 
This discrepancy in R e /r v between NFW and R 1 ' 4 law gets 
even worse if one adopts i? cu t = 4, which provides the best 
fits of the NFW surface mass density to the R 1 ^ 4 law: indeed, 
Figure || indicates _R C ~ 0.8 r v for the NFW model. 

In summary, the NFW surface mass density profile re- 
sembles an R 1 / 4 law in a fairly wide range of radii, but 1) 
one has to resort to an abnormally large effective radius, 
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Fig ure 12. The best fitting parameters of the Sersic law, eq. 
(|53|), as functions of concentration: 1/m (dashed line) and Re/r v 
(solid line). 

very close to the virial radius, and assume that the effective 
radius measures half the projected light (or mass) within 4 
times the virial radius, and 2) the fit is good in a consider- 
ably smaller range of radii than is observed in the nearby 
giant elliptical NGC 3379. 

The generalization of the i? 1//4 law into an R}' m law, 
first proposed by Sersic (1968), is known to fit the surface 
brightness profiles of elliptical galaxies within a much larger 
mass range than the de Vaucouleurs law (Caon, Capaccioli & 
D'Onofrio 1993). The surface brightness of the Sersic profile 
is 

I(R) = h exp{-b(™) [(R/R c ) 1/m - 1]} , (53) 

where b(m) is tabulated by Ciotti (1991), who gives the em- 
pirical relation b(m) ~ 2 m — 0.324, good to 0.1% relative 
accuracy. The de Vaucouleurs law is reproduced for m = 4, 
while m = 1 corresponds to an exponential law as in spiral 
disks. 

In Figure [H], we plot the NFW surface brightness 
I = Em/T, with Em given by equation ( fll| ) and c = 10, 
as a function of (R/R c ) 1/m for various values of the Sersic 
parameter m. We compare them to the Sersic profiles given 
by the straight dashed lines. The agreement is good for all 
values of m, within ranges of I/I e that increase with increas- 
ing m. Comparison of the plots for different m shows that 
the Sersic models with lower m generally agree better with 
the NFW surface brightness for smaller radii, while those 
with larger m are in better agreement at larger radii, closer 
to the virial radius. Overall, the NFW profile matches best 
the m = 3 Sersic law, over a factor of 10 3 in intensity (7.5 
magnitudes) . 

For a more quantitative comparison, we performed two- 
parameter fits of the Sersic models (B3J) to the projected 
NFW formula ([«]). The NFW profile was sampled in the 
range of 0.01 < R < 1 with a given c. The fitted Sersic pa- 
rameters 1/m and _R e /r v obtained for different c are shown 
in Figure uA Figure n3l compares the two projected profiles 



0.01 0.1 1 

R/r v 

Figure 13. Comparison of the projected NFW density profile, 
eq. ([H]), and the best fitting Sersic law, eq. ([^). 

for c = 10. The best-fit parameters of the Sersic model in 
this case are m = 3.07 and R c /r v = 0.55. 

While Caon et al. (1993) find similar ranges of agree- 
ment between observed profiles and Sersic laws, this range 
in intensity is still smaller than the range of 10 4 found for 
NGC 3379 by de Vaucouleurs & Capaccioli (1979). More- 
over, while Caon et al. (1993) find that the best fitting Sersic 
models for elliptical galaxies have indices spanning a wide 
range, from m = 2 for faint ellipticals to m = 10 for bright 
ellipticals, the Sersic laws that match the NFW models span 
a much smaller range, roughly m = 3±0.5 (2.71 < m < 3.41 
for 5 < c < 15, see Figure h2J). Moreover, the problem of very 
high values of R c /r v (0.46 < R c /r v < 0.81 for 5 < c < 15, 
see Figure ^), remains in the fits of Sersic profiles to pro- 
jected NFW models. 



4 DISCUSSION 

The main disadvantage of the NFW model is the logarithmic 
divergence of its mass (and luminosity for constant mass- 
to-light ratio). In contrast, the Jaffe (1983) and Hernquist 
(1990) models converge in mass, and their properties can be 
expressed in units of their asymptotic mass. For the NFW 
model, one is restricted to a mass at a physical radius such 
as the virial radius. This mass divergence also complicates 
the analysis of surface brightness profiles, which involve the 
effective radius where the aperture luminosity is half its 
asymptotic value. However, independently of the radial cut- 
off introduced to define the effective radius, the projected 
NFW density profile is consistent with constant mass-to- 
light ratio, given the observed Sersic profiles of elliptical 
galaxies, but only in a limited range of radii, with unusually 
high values of R c and in a smaller interval of Sersic shape 
parameters than observed. On the other hand, the Hern- 
quist (1990) model, whose density profile scales as r -4 at 
large radii, produces better fits to the R 1 ^ 4 law. 

The upper panel of Figure [Io| suggests that, for reason- 
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able effective radii, if indeed dark matter follows the NFW 
profile, the mass-to-light ratio, T, is not constant but in- 
creases with radius, not only in the outer regions, as is in- 
ferred from the commonly accepted picture of galaxies em- 
bedded in more spatially extended dark haloes, but also in 
the inner regions. This is at odds with the observed kine- 
matics of ellipticals that Bertola et al. (1993) inferred from 
observations of ionised and neutral gas around specific el- 
lipticals. Moreover, increasing T throughout the galaxy im- 
plies radial velocity anisotropy throughout elliptical galax- 
ies, whereas violent relaxation should cause isotropic cores.Q 
Thus it appears difficult to reconcile the photometry and 
kinematics of elliptical galaxies with NFW models. In a 
forthcoming paper (Mamon & Lokas, in preparation), we 
will omit the assumption of mass follows light in a more de- 
tailed assessment of the compatibility of the observations of 
elliptical galaxies with the NFW model. 

The results presented in this paper can be directly ap- 
plied to the analysis of the mass and light distribution in 
clusters of galaxies. A standard procedure to do it is to mea- 
sure the surface brightness and the light-of-sight velocity dis- 
persion and assuming some form of velocity distribution or 
mass-to-light ratio calculate the luminosity density and the 
velocity dispersion by solving the Abel integral equations 
( [fl| ) and (|ij) and the Jeans equation (Binney & Mamon 
1982, Tonry 1983, Solanes & Salvador-Sole 1990, Dejonghe 
& Merritt 1992). The results of this procedure are uncertain 
because it involves derivatives of observed quantities which 
are usually noisy. One also experiences a degeneracy because 
different models fit the data equally well (Merritt 1987). In- 
stead of solving the Abel equations one can also model the 
luminosity density and velocity dispersion with simple func- 
tions and fit their parameters so that they reproduce their 
projected counterparts (Carlberg et al. 1997). 

Our results are useful for the simpler approach of assum- 
ing realistic forms of the density distribution, velocity dis- 
tribution and mass-to-light ratio. Here we provide the tools 
for modelling the NFW density profile with different velocity 
distributions and constant mass-to-light ratio (T = const), 
and obtain exact predictions for the surface brightness and 
the line-of-sight as well as aperture velocity dispersion that 
can be directly compared to observations. 
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